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Physical properties of DvPm have been studied by magnetization, specific heat, and ultrasonic 
measurements. The magnetic entropy change and ultrasonic properties in intermediate phase 
11 indicate that the degeneracy of internal degrees of freedom is not fully lifted in spite of the 
formation of magnetic order. The ultrasonic attenuation and huge softening of C44 in phase 
II suggest the existence of electric-quadrupolar (orbital) fluctuations of 4/-electrons. These 
unusual properties originate from a geometrical quadrupolar frustration. 
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A geometrical frustration, that is, a geometrical prohi- 
bition against satisfying the lowest-energy state of each 
bipartite bond, leads to various unusual physical prop- 
erties, that often cause an extensive degeneracy in the 
ground state of a system and prevent ordering. A well- 
known example is the spin frustration in triangular anti- 
ferromagnets, such as CsCoCbj and CsNiChj, that causes 
various magnetic-ordering processes. 1-5 ' A renewed in- 
terest in geometrical frustration systems has resulted 
from the discoveries of new magnetic properties in cer- 
tain materials. Two examples of these materials are the 
pyrochlore oxides Dy 2 Ti 2 07 and Ho 2 Ti 2 Oy considered 
to be model systems of spin ice materials. 6 ' 7 - ) A spin 
ice model shows a strong geometrical frustration arising 
from the competition between a ferromagnetic nearest- 
neighbor interaction and an axial single-ion anisotropy. 
The spins in this model are disordered with local configu- 
rations obeying the ice rule in every tetrahedron, and the 
ground state is macroscopically degenerated with a resid- 
ual entropy. In fact, Dy 2 Ti 2 07 shows a residual entropy 
of ~ 1.68 J/mol K, which is numerically in agreement 
with Pauling's entropy for water ice, R(l/2) ln(3/2). 8 ) 
Another new noteworthy discovery is the formation of 
unusual short-range order named the magnetic random 
layer (MRL) lattice in RB 2 C (R = Dy, Ho), 9 ) that results 
in an unusual neutron scattering in the low-wave-numbcr 
region. The unconventional MRL state with large mag- 
netic fluctuations originates from the frustrated nature 
of the Shastry-Sutherland lattice (SSL) consisting of the 
rare-earth atoms. 10 ) 

In this letter, we focus on a frustration of the orbital 
(quadrupole) degrees of freedom for 4/ electron system, 
and propose that DyB4 is the first example of a geometri- 
cally quadrupolar (orbital) frustrated system. DyB4 crys- 
tallizes in the tetragonal ThLVtype structure belonging 
to space group D\ h — Pi/mbm. 11 ^ 13 ^ The boron atoms 
form a continuous three-dimensional network. The Dy 
atoms are located below and above the centers of the 
seven-membered boron rings in densely packed planes 
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of boron atoms. All Dy atoms occupy the equivalent 
sites having orthorhombic symmetry C 2l , (Fig. 4 (a)). 
The two-dimensional layer of the Dy atoms is formed by 
fused equilateral triangles and squares. It can be con- 
sidered that this arrangement of the Dy atoms is the 
SSL. Previous studies of the magnetic properties of DyB 4 
have shown that it orders antiferromagnetically. 14, 15 ) 
Two second-order magnetic phase transitions have been 
discovered at 20.4 K and 12.8 K. 16) A neutron diffrac- 
tion study of Dy n B 4 was performed at 4.2 K. 17 ) The 
deduced magnetic structure is a simple collinear antifer- 

romagnetic one with the sequence (+H ) for the four 

Dy atoms in the tetragonal unit cell. The magnetic mo- 
ments are oriented parallel to the [001]-axis. However, 
it seems that the data is insufficient for determining the 
detailed magnetic structure of DyB 4 , owing to the strong 
absorption of thermal neutrons by dysprosium. There are 
very few reports on the physical properties of DyB 4 and 
studies from the point of view of the frustration and the 
role of quadrupolar interactions are lacking. Here we in- 
vestigate the details of the physical properties of DyB 4 
by means of magnetic susceptibility, specific heat, and 
ultrasonic measurements. 

Single crystals of DyB 4 were prepared by the Czochral- 
ski method in argon atmosphere using a tetra-arc fur- 
nace. The magnetization of DyB 4 was measured using a 
SQUID magnetometer in the temperature range of 1.8- 
300 K. The magnetic susceptibility measurements were 
performed under zero-field-cooled (ZFC) and field-cooled 
(FC) conditions. Specific heat measurement was carried 
out by the heat-relaxation method down to 1.9 K. The 
elastic properties of DyB4 were investigated by means 
of ultrasonic measurements. Ultrasonic measurement is 
an effective technique for examining the quadrupolar ef- 
fect in an /-electron system with an orbitally degenerate 
ground state, because the quadrupole of the 4/ electron 
couples to the elastic strain. The elastic constant Cr of 
the sample was obtained by Cr = pv 2 , where v and p rep- 
resent the sound velocity and the mass density, respec- 
tively, v was measured by a phase comparator using dou- 
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Fig. 1. Magnetic susceptibility \ of DyB4. The inset shows the 
Curie- Weiss fit to the inverse susceptibility. 
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Fig. 2. Magnetic specific heat divided by temperature C mag /T 
(left axis) and corresponding entropy S (right axis) vs temper- 
ature for DyP>4. The horizontal broken lines indicate the values 
of entropy for one and two Kramers doublets. 



ble balanced mixers and p = 6.783 x 10 3 kg/m 3 was cal- 
culated using the lattice constants a = 7.0960(4) A and 
c = 4.0128(2) A of DyB 4 at 25 K. 18 ) 

Figure 1 shows the magnetic susceptibility x for a 
single crystal of DyB 4 . No difference between the FC 
and ZFC susceptibilities in any direction is observed 
throughout the whole temperature range; therefore, only 
ZFC data are shown for clarification. The coincidence 
of the ZFC and FC data suggests that no spontaneous 
magnetization is present in DyB 4 . The susceptibility for 
H ]| [001] shows a maximum at about 21.0 K, and a 
small discontinuous drop occurs at 20.3 K. In addition, a 
small shoulder is observed at 13.0 K. On the other hand, 
susceptibilities for H || [100] and H || [110] show no 
anomalies around 21 K, but exhibit cusplike maxima at 
13.0 K. These results suggest that the z and xy compo- 
nents of the magnetic moments order independently at 
T C i ~ 20.3 K and T C2 ~ 13.0 K, respectively. The results 
may be divided into three temperature regions, which are 
hereafter called phase I (T > T C i), II (T C2 < T < T C i) 
and III (T < Tq2) in the sequence of decreasing temper- 
ature. The magnetic structure in phase II is a partially 
dynamical one, that is, the components of the magnetic 
moments parallel to the [001]-axis are ordered while the 
perpendicular components remain disordered; then, be- 
low Tc2, all components are ordered. 

The inverse molar susceptibilities 01 DyB 4 (the 
inset of Fig. 1) show linear temperature dependences 
above 250 K for H || [001] and 35 K for H || [100] 
and H || [110]. The effective paramagnetic moments 
jUeff and the Weiss temperatures 9 P resulting from the 
Curie-Weiss fitting are as follows: p)^ 01 ' = 10.44 p B 
and 9l [wl] = -8.46 K; p^ m] = 10.90 p B and tff 001 = 

-33.40 K; p l ^ W] = 10.75 p B and 9 p lm = -31.42 K. 
Since the effective paramagnetic moments are close to 
that expected for a free Dy 3+ ion, i.e., 10.63 ps, we con- 
clude that the 4/ electrons in DyB 4 are well localized. 
The antifcrromagnetic 8 p are comparable to T"ci- These 
indicate that the spin frustration is small in DyB 4 . 

The specific heat of the single-crystalline sample of 
DyB 4 shows two A-likc anomalies that indicate exact 
transition temperatures of Tqi — 20.3 K and Tc2 = 
12.7 K. These transition temperatures are consistent 



with those reported previously 16 -' Since the peak of the 
specific heat at Tci is very sharp, and moreover, \ for 
H || [001] shows the discontinuity at Tci, the phase tran- 
sition at Tci seems to be a first-order one. In our present 
study, however, we cannot conclude whether this phase 
transition is a first-order or a second-order transition be- 
cause the discontinuity of \ at Tci is very small and the 
elastic constants are continuous at Tci (Fig. 3). Only 
the z components of the magnetic moments order at Tci; 
that is, the Z2-symmetry of the system breaks in phase II. 
Therefore, the phase transition at Tci may be a second- 
order one and may belong to the three-dimensional Ising 
universality class in which the peak of the specific heat 
at a phase transition is very sharp. The phase transi- 
tion at T32 is a typical second-order one. The tempera- 
ture dependence of the magnetic specific heat divided by 
temperature C mag /T is shown in Fig. 2. The estimated 
magnetic contribution C mag to the specific heat of DyB 4 
is obtained by subtracting the phononic part from the 
total specific heat of DyB 4 ; the phonon contribution is 
obtained from a measurement of the isostructural and 
nonmagnetic compound LuB 4 . We obtain the tempera- 
ture dependence of the magnetic entropy per Dy 3+ ion 
S by numerically integrating the data of C mag /T vs T. 
The magnetic entropy of approximately R In 2 and R In 4 
is released below Tc2 and Tci, respectively. The CEF 
ground-state multiplet J = 15/2 of Dy 3+ splits into eight 
E 1 / 2 (^5) doublets in the CEF potential of the site with 
Ci v symmetry in DyB 4 . The result indicates that the 
CEF ground state of DyB 4 is a pseudo-quartet consist- 
ing of two Kramers doublets with E 1 / 2 (^5) symmetry It 
is noted that the degeneracy due to the internal degrees 
of freedom is not fully lifted in spite of the formation of 
magnetic order in phase II. A similar entropy change is 
shown in the quadrupolar compound DyB2C2- 19 - ) In this 
case, the twofold Kramers degeneracy is conserved since 
time-reversal symmetry is not broken by the antifcrro- 
quadrupolar ordering. On the other hand, in a conven- 
tional magnet, degeneracy is lifted below the magnetic 
transition temperature since Zeeman splitting caused by 
an internal magnetic field is different in each Kramers 
doublet. Hence, phase II of DyB 4 is not the conventional 
magnetically ordered state. 
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Fig. 3. Temperature dependence of the elastic constants of DyB4. 
Transverse modes (Cn — Ci2)/2, C44, and C66 were measured 
at frequencies of 8 MHz. The insets show the low-temperature 
data of the transverse modes. The solid lines are fits to the data 
using C r (T) = C°(T - T°)/(T - 0). 



Elastic constants of DyB4 (Fig. 3) show clear anoma- 
lies associated with the successive transitions at Tci and 
Tc2- The transverse mode C44 associated with elastic 
strain e yz or e zx exhibits a large softening of 10.4% be- 
low 120 K down to Tci. On the other hand, the trans- 
verse mode (Cn — Cyi)j1 and Cqq show small softening 
of 0.8% and 0.6% between 70 K and Tci, respectively. 
The largest softening of C44 in phase I indicates that the 
dominant components of quadrupolc moments in DyB4 



are O 



yz 



JyJ z 



J z J y and O z 



J Z J X + J X J Z which 



couple with the elastic strains e yz and e zx , respectively. 
Furthermore, in phase I, the elastic softenings are de- 
scribed by the equation C T (T) = C$(T - Tg)/(T - Q), 
which is deduced from the Curie term of the quadrupolc- 
strain susceptibility. Here, 9 is proportional to the av- 
erage quadrupole-quadrupole interaction between the 
Dy 3+ ions. = 9 + Tj T , where E JT is the Jahn- 
Teller (4/ electron-lattice) coupling. 20 ) The character- 
istic temperatures obtained are as follows: 9 = 2.26 K 
and T° = 4.61 K for C 44 ; 6 = 3.95 K and T° = 4.19 K 
for C 66 ; 6 = 7.12 K and T° = 7.48 K for (C n - C 12 )/2. 
The positive values of 6 indicate that the predominant 
quadrupolar interactions in DyB4 are ferro-type. C44 
shows a kink at Tci and decreases abruptly with decreas- 
ing temperature from Tci to Tc2- It should be noted that 
C44 shows a marked softening of 40% and its echo signal 
shrinks considerably throughout the whole temperature 
range of magnetically ordered phase II. With decreasing 
temperature, the elastic constants change their behav- 
ior at T32 from softening to hardening, and then almost 
level off towards K. The echo signal is recovered and 
stabilized below Tc2- 

The huge elastic softenings and the ultrasonic absorp- 



tion of C44 in phase II are different from the critical 
behaviors observed in the vicinity of a structural phase 
transition temperature induced by either ferroquadrupo- 
lar ordering (e.g., DyB@, HoBg 21 ') or the cooperative 
Jahn- Teller effect (e.g., DyV0 4 22) ). On the other hand, 
it has been known that large elastic softenings and ultra- 
sonic absorptions in an ordered state arise from magnetic 
and structural domain-wall stress effects. 23 ) However, in 
DyB 4 , no symmetry reductions of the crystal and mag- 
netic structure at the transition from phase II to phase 
III were observed in the neutron diffraction study. 18 ) 
Therefore, the elastic instability observed in only phase II 
of DyB4 also cannot be explained by domain- wall stress 
effects. Hence, we conclude that the attenuation of the 
echo signal of elastic constants in phase II of DyB4 is 
caused by quadrupolar fluctuations and that the orbital 
degeneracy is still conserved in phase II. This is consis- 
tent with the results of the magnetic entropy change. On 
the contrary, the elastic behaviors in phase III of DyB 4 
indicate that the degeneracy of the orbital degrees of 
freedom is lifted and the quadrupolar fluctuations disap- 
pear. Consequently, quadrupolar ordering and antiferro- 
magnetic ordering may coexist below Tc2- 

Hereafter, we focus on the quadrupolar instability in 
phase II of DyB4, and we hypothesize that the quadrupo- 
lar fluctuations in DyB4 can be caused by a frustration in 
the 4/ orbital (quadrupolar) system of DyB4. In the fol- 
lowing, we consider two possible origins of the quadrupo- 
lar frustration in DyB4. The first is based on a competi- 
tion between the quadrupole-quadrupole interactions on 
the SSL formed by Dy 3+ ions in the xy-layer. An ele- 
mentary unit model of the SSL for the nearest-neighbor 
interactions is shown in Fig. 4 (c). If the coupling K 
is of antifcrro-typc, the quadrupolar system is frustrated 
regardless of the sign of the coupling K' , as in a spin sys- 
tem. The second is that the uniaxial anisotropy effects of 
the 4/ orbitals compete with the ferroquadrupolar inter- 
actions. The single ion anisotropy of the 4/ orbital prob- 
ably reflects the local environment of the rare-earth site 
through 4/ electron-lattice coupling and/or CEF effects. 
The local symmetry C^ v of the Dy site is characterized by 
a twofold axis (z'-axis in Fig. 4 (a)) lying in the xy-plane. 
Because the site symmetry of the Dy site is uniaxial, it 
is reasonable to assume that the quadrupolar anisotropy 
is also uniaxial. The most important point to note is the 
geometrical arrangement of the twofold axis at the Dy 
site. In Figs. 4 (b) and 4 (c), the twofold axes at the Dy 
sites are represented by arrows. The two twofold axes of 
a Dy-pair coupled by the X-bond arc collincar and the 
directions are antiparallel to each other. In contrast, the 
arrows between the pairs are orthogonal to each other. 
If only the single-ion anisotropy of Dy 3+ is considered, 
4/ orbitals on the neighboring sites are orthogonal to 
each other. Therefore, the single-ion anisotropy effects of 
the 4/ orbital in DyB4 compete with the ferroquadrupo- 
lar interactions. In addition, the present results that O 
are comparable to Ejt suggest that this competition is 
strong in DyB4. Accordingly, it can be expected that the 
quadrupolar system of DyB4 is geometrically frustrated 
and shows a quadrupolar instability. It can be said that 
this quadrupolar frustration model is similar to the spin- 
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Fig. 4. (a) Local structure of Dy 3+ ion in DyB<i. The Dy 3 + ion 
is indicated by the gray circle and the boron atoms by the open 
circles, (b) Lattice structure of the Dy 3+ ions of DyB<i. The 
bold solid lines are the inter square-lattice bonds K, and the 
bold broken lines are the intra square-lattice bonds K' . The thin 
solid lines represent the boron frame. The twofold axes of the 
local symmetry at the Dy sites are represented by the arrows, 
(c) Interactions of the elementary unit of the SSL. 



ice model which shows that strong geometrical frustra- 
tion can arise in presence of the dominant ferromagnetic 
interactions and local Ising anisotropy. 6, 7 > We obtained 
the result that predominant quadrupolar interactions in 
DyB4 are of ferro-type. However, these results can be ob- 
tained if the coupling K' and the interlayer interactions 
are ferro-type because 9 is only a mean value of the in- 
teraction. At present, the information is insufficient to 
conclude which origin of the frustration is more credible. 

In the following, we compare DyB4 and RB 2 C (R = 
Dy and Ho) in which the rare-earth ions also form the 
SSL. DyB 2 C and HoB 2 C show a three-dimensional mag- 
netic order but MRL states still remain below T c . 10 ) The 
magnetic ordering temperatures of RB 2 C (R = Dy, Ho, 
and Er) do not follow the de Gennes relation. This fact 
suggests that the ordering temperatures of DyB 2 C and 
HoB 2 C are suppressed. The unconventional fluctuating 
ground states and the suppressions of the ordering tem- 
peratures of DyB 2 C and HoB 2 C are due to spin frustra- 
tion effects. In DyB4, however, the large magnetic fluc- 
tuations such as those observed in DyB 2 C and HoB 2 C 
were not recognized from the results of powder neutron 
diffraction measurement 18 -' ; this is consistent with the 
results of the inverse magnetic susceptibilities. Conse- 
quently, it is thought that the influence of spin frustration 
on the magnetic transition temperature is very small in 
this compound. It seems that the quadrupolar system of 
DyB4 is strongly influenced by a geometrical quadrupo- 
lar frustration and that the quadrupolar ordering is sup- 
pressed. 

Most of the compounds with spin frustration show 
long-range order in the lowest-temperature phase. How- 
ever, a spin frustration system increases the degree of 
order through successive phase transitions accompanied 
by high-entropy intermediate states with spin fluctua- 
tion, because the energy gain of the spin-spin correla- 
tion is suppressed by the spin frustration. For example, 
a partial disorder and an order of limited components of 
spins are generally observed in triangular lattice antifcr- 
romagnets such as CsCoCl 3 and CsNiCl 3 . 1_5) CsNiCl 3 
exhibits successive phase transitions similar to those of 



DyB 4 , where the parallel and perpendicular components 
to the c-axis of Ni 2+ moments order at Tni = 4.8 K and 
Tn2 = 4.4 K, respectively. In this material, the large dif- 
ference between the Weiss temperature 8 P = — 69 K and 
Tni is attributed to the spin frustration. 5 -* On the other 
hand, in DyB 4 , the successive phase transitions charac- 
terized by a magnetic ordering of the z and xy com- 
ponents at independent temperatures are not caused by 
spin frustration. The quadrupolc moments and the xy 
components of the magnetic moments order at Tq 2 si- 
multaneously. The quadrupolar interactions most likely 
play a dominant role in the phase transition at Xc2, and 
the ordering of the xy components of the magnetic mo- 
ments coupled with the quadrupole ones can be induced 
by the ordering of the quadrupole moments. At present, 
however, it is unclear why the phase transition at Tc 2 is 
not a first-order but a second-order transition. Accord- 
ingly, we conclude that the successive phase transitions 
in DyB 4 are caused by the coexistence of quadrupolar 
frustration and the quadrupolar entropy effect. 

In summary, we observed successive phase transition at 
T C i = 20.3 K and T C2 = 12.7 K in DyB 4 . The successive 
phase transitions are characterized by a magnetic order- 
ing of the z and xy components at independent temper- 
atures. The magnetic entropy change indicates that the 
degeneracy of the internal degrees of freedom in phase 
II is still conserved in spite of the formation of magnetic 
order. The huge softening of C44 and the attenuation in 
phase II suggest the existence of quadrupolar fluctua- 
tions. These unusual properties originate from a geomet- 
rical quadrupolar frustration. 

We are grateful to T. Sakakibara, J. Custers, and 
N. Todoroki for their helpful discussions. 
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